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Abstract. Monte Carlo techniques have been widely employed in statistical physics as well as in quantum theory in the 
Lagrangian formulation. However, in the conventional approach, it is extremely difficult to compute the excited states. 
Here we present a different algorithm: the Monte Carlo Hamiltonian method, designed to overcome the difficulties of 
the conventional approach. As a new example, application to the Klein-Gordon field theory is shown. 



INTRODUCTION 

There are two standard formulations in quantum the- 
ory: Hamihonian and Lagrangian. A comparison of the 
conventional approaches is given in Tab. 1 . 

Monte Carlo (MC) method with importance sampling 
is an excellent non-perturbative technique to calculate 
path integrals in quantum theory. In the last two decades, 
and it has successfully been applied to Lagrangian lattice 
gauge theory (1, 2, 3). In the standard Lagrangian MC 
method, however, it is extremely difficult to compute the 
spectrum and wave function beyond the ground state. On 
the other hand, the standard Hamiltonian formulation is 
capable of doing it. 

Recently, we proposed an algorithm to construct an ef- 
fective Hamiltonian from Lagrangian MC simulations in 
Ref. (4). We called it the MC Hamiltonian method. The 
advantage, in comparison with the standard Lagrangian 
MC approach, is that one can obtain the spectrum and 
wave functions beyond the ground state. It also allows to 
do thermodynamics. In this paper, we briefly review what 
we have done, and present some new results. 



ALGORITHM 
Effective Hamiltonian 

Let us review briefly the basic ideas of our approach. 
The (imaginary time) transition amplitude between an 
initial state at position x,, and time f,, and final state at 
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Table 1. Comparison of the conventional methods. 



Formulation 


Hamiltonian 


Lagrangian 


Approach 


Schrodinger Eq. 

H\Er, >= E„\E„ > 


Path Integral 

<0> = 

Jd[x]0[x]sxp{-Slx]/ri) 






/[rf.v]exp(-S[x]//!) 


Algorithm 


Series expansion, 
variational approx., 
Runge-Kutta ... 


MC simulation 
with importance 
sampling 


Advantage 


Both ground state, 
& excited states 
can be computed. 


It generates 
the important 
configurations 


Problem 


Difficult for 
many body systems. 


Difficult for 
excited states. 



Xf and tf is related to the Hamiltonian H by 

Mfi ^<xf,tf\xi,ti >=<x/|e-^('/-'')/'^|xi > 

oo 

= ^ < Xf \E„ > e-""^^" < E„ \xi > , (1) 

,1=1 

where T — tf — ti. According to Feynman's path integral 
formulation of quantum mechanics (Q.M.), the transition 
amplitude is also related to the path integral: 

Mfi^ l[dx]exp{~S[x]/h)f,^:;f. (2) 

The starting point of our method, as described in detail 
in Ref. (4) is to construct an effective Hamiltonian Hgff 
(finite N xN matiix) by 

Mfi^<Xf\e-"'ff^/^'\xi> 



= I < ^f\E'J^ > e-"""^'^ < Kff\xi > . (3) 

n=l 

The eigenvalues E^^^ and wave function > can be 
obtained, by diagonalizing M using a unitary transforma- 
tion 



M = U^DU, 



(4) 



where D = diag{e^'^i^^^^, ...,e^^N ^A). Once the spec- 
trum and wave functions are available, all physical infor- 
mation can also be obtained. 

Since the theory described by H, which basis in 
Hilbert space is infinite, is now approximated by a theory 
described by a finite matrix H^ff, which basis is finite, the 
physics of H and He/f might be quite different at high en- 
ergy. Therefore we expect that we can only reproduce the 
low energy physics of the system. This is good enough 
for our purpose. In Refs. (4, 5, 6, 7), we investigated 
many 1-D, 2-D and 3-D Q.M. models (Tab. 2) using this 
MC Hamiltonian algorithm. We computed the spectrum, 
wave functions and some thermodynamical observables. 
The results are in very good agreement with those from 
analytical and/or Runge-Kutta methods. 

Table 2. Q.M. systems, investigated by the MC Hamiltonian 
method using the regular basis. 



System 


Potential 


Q.M. in 1-D 


V{x)=0 
V{x) = imroV 
V{x) = -Vosech2(x) 

V{x) = ^\x\ 

^ Fx, x>0 


Q.M. in 2-D 


V(x,y) = ^mOi^x^ + jmG?y^ 
V{x,y) — jma?x^ + jinap-y^ + hey 


Q.M. in 3-D 


V{x,y,z) = jinoy^x^ + jm(£p-y^ + jiMO^z^ 



Basis in Hilbert Space 

To get the correct scale for the spectrum, the position 
state \xn) (Bargman states or box states) at or tf should 
be properly normalized. We denote a normalized basis of 
Hilbert states as l^n), n — l,...,iV. In position space, it 
can be expressed as 



e„(x) = { 



l/^/AXn, x€[Xn,Xn+\] 
0, X ^ [Xn,X„+i] 



(5) 



where Ajc„ = x„+ 1 — x„ . 

The simplest choice is a basis with Ax„ = const., 
which is called the " regular basis". In Refs. (4, 5, 6, 7), 
the regular basis is used. For many body system or quan- 
tum field theory, the regular basis will encounter problem. 
For example, in a system with a 1-D chain of oscillators 
(see later), if the number of oscillators is 30, the mini- 
mum non-trivial regular basis is A'^ = 2-'*' = 1073741824, 
which is prohibitively large for numerical calculations. 

Guided by the idea of important sampling, in Refs. 
(8, 9), we proposed to select a basis from the Boltzmann 
weight proportional to the transition amplitude between 
X- = at /■ = and x'j- = x„ at some t'j^. In a free parti- 
cle or harmonic oscillator case, the distribution is just a 
Gaussian 



Pbasis [Xn] 



1 



exp(-;^4) 
27ta ^ 2a2^ 



(6) 



where a = ^ Ht'r/m for the free case and a 



'/zsinh(cofy-)/(mCo) for the harmonic oscillator. We call 
such a basis the "stochastic basis". 

Matrix elements 

As explained above, the calculation of the transition 
matrix elements is an essential ingredient of our method. 
The matrix element in the normalized basis is related to 

<x„i,tf\xn,ti > by 



< e„/,r|e„,0 > 

r^'dx'f 



W V'Ax„'Ax„ < Xni,tf\x„,ti >, (7) 

where < x„',?/|x„,?, > can be calculated using MC simu- 
lations as follows. 

(a) Discretize the continuous time. 

(b) Generate free configurations [x] between t G {ti,tf) 
obeying the Boltzmann distribution 

exp(-5o[x]//z) 



Po[x] = 



/[Jx]exp(-5o[x]//z) ' 



where So = Jdt m£- /2. 
(c) Measure 

<0v>= j [dx]cxp(-Jdt vix)/n)\l"^;lfPo[x] 

The path integral in Eq. (7) is then 

, ^ / m 

<X„>,tf\Xn,ti> = <Ov>X\ll^^ 

X exp[- — (x„/-x„) 



(8) 



(9) 



(10) 



QUANTUM FIELD THEORY 

The main purpose of the algorithm is to study many 
body systems and quantum field theory beyond the 
ground state. As an example, we consider a chain of Nose 
coupled oscillators in 1 spatial dimension. Its Hamilto- 
nian is given in Ref. (10) as 



(11) 



7=1 



where pj and qj are the momentum and displacement of 
the j-th oscillator respectively. This model is equivalent 
to the Klein-Gordon field theory on a (l+l)-dimensional 
lattice. 

The spectrum of the system is analylically known: 

1, 



a2(2sinfe/2)2 + i2g, 



(12) 



where ni,...,nNg^^ = 0, 1, . . ., and k = 2%1/Nosc with / an 
integer between —Nose/ '2 and Nose/2. 

We generate a stochastic basis according to Eq. (6) 
with N=1000 configurations [qi, ■■■,qNosc] for the initial 
and final states iorNose = 9 oscillators. For the adjustable 
parameter o, we use the distribution for the uncoupled 
oscillators at /y^ = // for simplicity. (Of course, one should 
study systematically the dependence of the results on a). 
Tab. 3 compares the spectrum from the MC Hamiltonian 
with the analytical results for the first 20 states with Q, = 
1,0.0 = 2, m= I, H= I, and T = 2. They agree very well 
with the exact ones. 



SUMMARY 

In this paper, we have tested the MC Hamiltonian 
method with a stochastic basis in a many body Q.M. sys- 
tem with a chain of coupled oscillators: the Klein-Gordon 
field theory on a (l+l)-dimensional lattice. The results 
are very encouraging. We believe that the appUcation 
of the algorithm to more complicated body systems and 
quantum field theory will be very interesting. 
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Table 3. Comparison of the spectrum of the chain 
of 9 coupled oscillators, between the MC Hamil- 
tonian method with a stochastic basis and the an- 
alytic ones. 
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1 ^ (\AA'X'\ 1 <Q0^/I'7 
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1 J.UD /<5UjcSoy4cS4 
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i:>.zyyyD / :)4iz4z 


\'x n 1 *^ni QQ^'Jcn 

1 j.jziDUiyyjjou 


o 


1^ 1A<AQr\i^1Q^C\A 


ij.jziouiyyjjou 


1 


1^ 0<1 1'S^QT 

ij.jjziyjijjos/ 


1 1 <QOQ 1 1 'TOI '7'2'J 

1 j.joyoi 1 /yi /j3 


8 


13.585794986361 


13.589811791733 


9 


13.680136748933 


13.751084748745 


10 


13.744919087477 


13.751084748745 


11 


14.984737011385 


14.944060480668 


12 


15.012353803145 


15.057803869484 


13 


15.057295761044 


15.057803869484 


14 


15.108904652020 


15.171547258300 


15 


15.125356713561 


15.171547258300 


16 


15.187413290039 


15.171547258300 


17 


15.308536490102 


15.321601993380 


18 


15.396255686587 


15.321601993380 


19 


15.420708031412 


15.435345382196 


20 


15.432823810789 


15.435345382196 
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